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Abstract 

We determine the most general solution of the five-dimensional vacuum Einstein equation, 
allowing for a cosmological constant, with (i) a Weyl tensor that is type II or more special in 
the classification of Coley et al, (ii) a non-degenerate “optical matrix” encoding the expansion, 
rotation and shear of the aligned null direction. The solution is specified by three parameters. 
It is locally isometric to the 5d Kerr-de Sitter solution, or related to this solution by analytic 
continuation or taking a limit. This is in contrast with four dimensions, where there exist infinitely 
many solutions with properties (i) and (ii). 


1 Introduction 

There is a long tradition of classifying solutions of the four-dimensional Einstein equation according 
to the algebraic type of the Weyl tensor. If the Weyl tensor is everywhere type II, or more special, 
then the solution is called algebraically special. In an algebraically special spacetime, the Einstein 
equation simplifies considerably and one can determine the explicit dependence of the metric on 
one of the coordinates. The Einstein equation then reduces to PDEs in 3 dimensions. The general 
solution of these PDEs is not known. However, it is clear from various special cases that the general 
solution involves arbitrary functions [T]. 

If one makes the stronger assumption that the Weyl tensor is (everywhere) of type D then much 
more progress can be made. The vacuum Einstein equation can be solved explicitly [2]. There are 
several families of solutions, each specified by a few parameters. These solutions include the Kerr 
solution. Performing this classification led to the discovery of a new vacuum solution: the spinning 
C-metric, describing a pair of rotating black holes being accelerated by cosmic strings. 

In this paper we will consider higher-dimensional solutions of the vacuum Einstein equation, 
allowing for a cosmological constant: 

Rab — ^9 ah- (1) 

The algebraic classification of the Weyl tensor has been extended to d spacetime dimensions by 
Coley et al. [3]. A solution is type II or more special in this classification if it admits a multiple 
Weyl aligned null direction (multiple WAND) i. The condition for £ to be a multiple WAND can be 
written [3] 

^''^[eCa]b[cd^f] = 0 , ( 2 ) 

where Cabcd is the Weyl tensor. In four dimensions, this is the same as the condition for i to be 
a repeated principal null direction. The Myers-Perry black hole solution [5] (“higher dimensional 
Kerr”) is known to have a Weyl tensor of type D [6l [TJ E] , which means that it admits two distinct 
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multiple WANDs. The Myers-Perry solution has been generalized to include a cosmological constant 
in five [9] and higher m dimensions. For simplicity, we will refer to these as “Kerr-de Sitter” 
solutions (for any value of A). These also have a Weyl tensor of type D [7]. 

By analogy with the four-dimensional case, one might expect there to exist a small number of 
families of solutions of ([T]) with a Weyl tensor of type D. Surprisingly, we will establish a stronger 
result in five dimensions: subject to one extra assumption, the Kerr-de Sitter solution is essentially 
the only solution with a Weyl tensor of type II. 

To explain the extra assumption, we note that any solution of ([T]) admitting a multiple WAND 
must also admit a geodesic multiple WAND [H] hence there is no loss of generality in assuming i to 
be geodesic. Recall that the expansion, rotation and shear of the null geodesic congruence tangent 
to i are defined as the trace, antisymmetric part, and traceless symmetric part of the (d — 2) x {d — 2) 
“optical matrix” 

Pij = (3) 

where rrii are a set of (d — 2) orthonormal spacelike vectors orthogonal to i. Our assumption is that 
Pij is non-degeneratelll We will prove: 

Theorem. Let {M,g) be a solution of the five-dimensional vacuum Einstein equation ([T]). Assume 
that {M,g) admits a geodesic multiple WAND i for which the 3x3 matrix pij is non-degenerate. 
Then one can define an affine parameter r along the null geodesics tangent to I such that the 
eigenvalues of pij are 1/r and l/(r ± ix) for some real function y constant along each geodesic. 
Furthermore: 


1. If X 7^ 0 and dx 7^ 0 then one can define local coordinates {u,r,x,x,y) such that the metric is 


ds^ = —2(dn -|- x^dy) 


dr -F H{r, x)(du -|- x^'dy) - Eq ( dx - ^dy ) -F Pdy 




where 


+r^X^ ( dx - ^dy ) + {r^ + x^) 


A 


( dx^ 


\Pix) 


+ P{x)dy^ 


H{r,x) = Ao - -(r -x ) - 


Mo 


2(r2 -|- x^) ’ 

Ax^ 


P{x) = Co-^-2AoX^- , 

4 

{Aq, pq,Co, Eq) are arbitrary real constants, and P{x) > 0- 


(4) 

(5) 

( 6 ) 


2. If X / 0 and dx = 0 then one can define local coordinates {u, r, x, y^ ,y^) such that the metric 


IS 


ds^ = —2(dtt -|- 2xA) 


dr -|- H{r){du + 2xA) — Eq ( dx -|- 


+r‘^X^ ^dx -b + X^)Kii{y)dy°'dy^, 


X^ 


A 


(7) 


where hap{y) (1 < a, /? < 2) is the metric on a 2d Riemannian manifold of constant curvature, 
A = Aa{y)dy°‘ is a 1-form such that dA is a volume form for this 2d manifold. 


H(r) = ^ - -(r^ -b x^) - 

2x4 8^ 2(r2+x2) 


Mo 


( 8 ) 


^Note that this does not depend on how the rrii are chosen. 
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and (x,^ 0 )-£'o) are arbitrary real constants (with x 7 ^ 0). The Ricci scalar of hap is 


= ( 9 ) 

X 

3. If X = 0 then one can define local coordinates {u,r,y^ ,y‘^,y^) such that the metric is 

~ ~ ^ ~ “ 2dudr + r‘^hij{y)dy^dy^, (10) 

where k G {—1,0,1}, fj-o is a constant, and hij{y) (1 < i,j < 3) is the metric on a 3d 
Riemannian space of constant curvature with Ricci scalar 6k. 


We will now make some remarks on the above theorem. 

For any of the above metrics d/dr is a geodesic multiple WAND. The Weyl tensor vanishes if, 
and only if, fiQ = 0. For yo 0 each solution is type D, i.e., there exists a second multiple WAND. 
The metric dH), has a scaling symmetry: for A 7 ^ 0 one can perform a coordinate transformation 


u . , 

u = — r = \r 
A 


X = Ax' 

and the metric in the primed coordinates takes the same form as 


X' y' 


( 11 ) 

but with rescaled constants: 


A' 

^ 0 - 


u' 


r' 


F ' - -^0 

-V' 


( 12 ) 


This shows that the solution is really a 3-parameter family^ 

We will show (in section HD that the solution @ is locally isometric to the 5d Kerr-de Sitter 
solution [9] with two unequal rotation parameters o The coordinates of (HD are closely related to 
the coordinates for the Kerr-de Sitter solution defined in Ref. m- Only when the parameters 
{Aq, fiQ,Co, Eq) lie within a certain set does the metric (HD describe a regular black hole solution. 
Values of the parameters outside this set can result in local metrics such as the “Kaluza-Klein 
bubble” spacetime of Ref. m which was shown to be algebraically special in Ref. |14jFI 

The metric m has a scaling symmetry analogous to m- Hence this solution is really a 2- 
parameter family. If > 0, one can perform a coordinate transformation to show that the solution 
is locally isometric to the Kerr-de Sitter solution with two equal, non-zero, rotation parameters. 
Solutions with < 0 can be regarded as analytically continued versions of the Kerr-de Sitter 
solution. 

The metric (jlOp is a generalized Schwarzschild metric written in outgoing Eddington-Finkelstein 
coordinates. Of course, this corresponds to the Kerr-de Sitter metric with vanishing rotation pa¬ 
rameters. 

Note that, in four dimensions, there are many solutions satisfying the assumptions of the above 
theorem. In 4d, the Goldberg-Sachs theorem implies that i is shear-free, which implies that pij is 
degenerate if, and only if, pij = 0 (this defines the Kundt family of solutions). So any algebraically 
special solution with pij 7 ^ 0 satisfies the assumptions of the theorem. As noted above, there is no 
simple explicit form for such solutions, and such solutions are known to involve free functions. 

^We choose not to use this symmetry to eliminate one of the constants because this requires consideration of various 
special cases e.g. if So yf 0 then one can rescale to set So = 1 but then the case So = 0 needs to be discussed separately. 

^Except in the special case A = ^0 = So = 0, Co = /ro, which can only be obtained as a limit of the Kerr-de Sitter 
solution. This is the only non-trivial case for which the second multiple WAND has degenerate (in fact vanishing) 
optical matrix. 

■^This spacetime was originally obtained by analytic continuation of the ’’Boyer-Lindquist” coordinates of the Myers- 
Perry solution. No analytic continuation of the coordinates is required to obtain it in the coordinates of 
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The non-degeneracy assumption on pij cannot be eliminated from our theorem. To see this, note 
that one can take the product of a 4d Ricci flat algebraically special solution with a flat direction 
to obtain a 5d Ricci flat solution admitting a geodesic multiple WAND. This solution will have 
degenerate pij because i does not expand along the flat direction. Obviously there are as many such 
solutions as there are 4d algebraically special solutions. For A 7 ^ 0 one can take a warped product 
to reach the same conclusion. 

Our result can be viewed as a new kind of uniqueness theorem for the 5d Kerr-de Sitter solution. 
It should be contrasted with the usual uniqueness theorem nais] for the 5d Myers-Perry black 
hole, which assumes the existence of a M x isometry group, asymptotic flatness, a regular 

horizon of spherical topology, and no topology outside the horizon. Our result assumes nothing 
about isometries or global structure and allows for a cosmological constant. 

The local nature of our result is similar to the uniqueness result for spacetimes with certain 
“hidden symmetries”, i.e., symmetries associated to Killing tensors rather than Killing vectors. 
Refs, [mils] proved that the Kerr-de Sitter solution (generalized to allow for a NUT charge [T^ ) 
is the most general d-dimensional solution of ([T]) admitting a “principal conformal Killing-Yano 2 - 
form”. Note that this result applies even for d = 4, where it yields a subset of the type D solutions. 
This is in contrast with our result, for which there exist infinitely many d = 4 solutions satisfying 
the assumptions of the theorem. 

Our result is also reminscent of the theorem that asserts that the Kerr solution is the unique 
stationary solution of the 4d vacuum Einstein equation with vanishing Mars-Simon tensor |19j . 
However our theorem does not assume stationarity. 

There have been several hints that higher-dimensional algebraically special solutions with non¬ 
degenerate Pij might be more rigid than their 4d counterparts. First, Ref. [20] considered higher 
dimensional Robinson-Trautman solutions, dehned by the existence of a null geodesic i with 0 7 ^ 
Pij oc 5ij (such I must be a multiple WAND). It was found that these solutions are considerably 
simpler for d > 4 than for d = 4. For d = 5, the only such solution is the generalized Schwarzschild 
metric m- 

Second, Ref. m investigated the possible existence of families of algebraically special solutions 
that contain the Schwarzschild solution. A solution “close” to the Schwarzschild solution in such a 
family would have non-degenerate pij. The approach of Ref. |21] was to consider linear perturbations 
of the Schwarzschild solution that preserve the algebraically special property. For d = 4 there are 
infinite families of such perturbations. But for d > 4 it was found that the only such perturbations 
that are regular on the orbits of spherical symmetry are perturbations corresponding to a linearisation 
of the Myers-Perry solution around the Schwarzschild solution. 

Third, Ref. [22| studied algebraically special solutions in d > 4 dimensions that (i) have non¬ 
degenerate Pij and (ii) are asymptotically flat. It was assumed that the curvature components can 
be expanded in inverse powers of an affine parameter r along the null geodesics tangent to It was 
found that such solutions are non-radiative, in contrast with the d = 4 case. 

Ref. [131 showed that 5d algebraically special solutions can be classified according to the rank of 
Pij. Our theorem determines all solutions for which pij has rank 3. Rank 0 dehnes the Kundt class 
of solutions which was studied in Refs [23l|23|. For rank 2 or rank 1, all solutions for which i is 
hyper surf ace-orthogonal {p[ij] = 0) were determined in Ref. |25j . We intend to return to the general 
(non-hypersurface-orthogonal) rank 2 and rank 1 cases in future work. 

We end this introduction with an outline of the proof of our theorem. The starting point for the 
proof is the recent demonstration |14] that the optical matrix of a geodesic multiple WAND in 5d 
can be brought to a certain canonical form by an appropriate choice of the basis vectors mj. This is 
the 5d analogue of the “shearfree” property that holds in 4d because of the Goldberg-Sachs theorem. 
As noted in Ref. mi, case 3 of the theorem follows immediately from combining this canonical form 
with the results of Ref. |2nj . 

For non-degenerate pij, the canonical form involves two unknown functions. We show how the 
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evolution equation for pij can be integrated to determine the dependence of these functions on an 
affine parameter r along the geodesics tangent to This determines the form of the eigenvalues of 
Pij as stated in the theorem. 

Next we complete {^, ruj} to a null basis {£, n,mj} where n is null and orthogonal to m,. After 
exploiting a residual freedom in the choice of ruj, we show how the “Newman-Penrose” and Bianchi 
equations can be integrated to determine the r-dependence of the basis vectors and hence the r- 
dependence of the metric. The r-dependence of the connection and curvature components is also 
fully determined. This calculation reveals that the Weyl tensor is necessarily of type D. 

The vanishing of certain connection components enables us to introduce local coordinates in a 
canonical way. After expressing our basis vectors in terms of these coordinates and using the results 
obtained previously we obtain a set of equations that can be integrated. At this stage, it becomes 
convenient to divide the analysis into two cases depending on whether y is constant or not. If dy 7 ^ 0 
then we use y as a coordinate and show that residual coordinate freedom can be exploited to make 
the solution independent of three of the remaining coordinates. Finally we solve for the dependence 
on y to obtain the solution (jl]). If dy = 0 then a similar procedure leads to the metric Q. 

This paper is organized as follows: section [2] contains the first part of the proof of the theorem 
in which we determine the connection and curvature components in a null basis. In section [3l we 
introduce coordinates and complete the proof of the theorem. Section |4] demonstrates how the 
metrics © and ([7]) are related to the Kerr-de Sitter solution. 

Notation 

We will perform most calculations in a null basis. Refs. [26tl27j developed a higher-dimensional ana¬ 
logue of the Newman-Penrose formalism used for calculations in such a basis. This was repackaged 
into a higher-dimensional analogue of the Geroch-Held-Penrose (GHP) formalism in Ref. [28]. We 
will follow the notation of Ref. [28] for the connection components and Weyl tensor components. In 
particular, we refer the reader to eqns. NP1-NP4, B1-B8 and G1-G3 of Ref. |28|, which lists all the 
Newman-Penrose and Bianchi equations satisfied by the connection and curvature components, as 
well as equations for the commutator of derivatives. 


2 Integration of GHP equations 

2.1 Canonical form for pij 

Gonsider an Einstein spacetime, i.e. a solution of ©, admitting a multiple WAND i. Introduce a 
null basis {i,n,mi}, i = 2,3,4. The multiple WAND condition is equivalent to the vanishing of the 
Weyl components of boost weights +2 and -|-1: 

^ij = 0, Tjjfc = 0. (13) 

Without loss of generality, i can be assumed to be geodesic m, i-e. 

Ki = 0. (14) 

Ref. [14] showed that the spatial basis vectors (i = 2,3,4) can be chosen so that the optical 
matrix pij of i takes one of the following forms: 

/ 1 a 0 \ / 1 a 0\ /la 0 \ 

b I —a 10 I 61 —a 1 0 I 61 —a —a^ 0 j (15) 

\0 D l + a^ ) \0 00 / \0 0 0 / 
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for functions a, b. For 6/0, these matrices have rank 3,2,1 respectively. Our assumption is that pij 
is non-degenerate so 

/ 1 o 0 \ 

Pij = b i -a 1 0 I , 6/0. (16) 

\ 0 0 l + a^ J 

We know all the GHP scalars with positive boost weight, both from the connection and the curvature: 

Hi, and we know the structure of pij, which has boost weight -|-1. Using the Newman- 
Penrose equations and Bianchi identities, we can determine the GHP scalars of boost weight negative 
and zero by systematically examining these equations from higher to lower boost weight. We now 
proceed to indicate the steps involved in this calculation. 

2.2 Choice of basis 

The form of the optical matrix above does not fix the basis uniquely because this form is preserved 
by null rotations about i and spins in the 2-3 directions. We can use this freedom to make some of 
the GHP scalars vanish. Consider a null rotation about i with parameters Zi |28] 

(17) 

(18) 
then gives 

(19) 


1 2 

i i, n n + ZiiTii — 2 ^ mi rrii — zi^, 

where z^ = ZiZi. This leaves pij unchanged (as Ki = 0) but Tj changes according to 

Ti^ Ti + PijZj. 

Since pij is non-degenerate, we can choose Zi to set r* = 0. Equation NP2 of Ref. 
= 0 [28]. In summary, we choose our basis so that: 

n = r'i= 0. 

It is convenient to combine the spatial vectors m 2 , m 3 into complex null vectors: 


m5 = 


m 2 -I- zms 


m5 = 


777-2 ~ ”^^3 


V2 


In this frame, we hav^ 



( 20 ) 


( 21 ) 


while the optical matrix is written as 



( 22 ) 


Now consider a spin in the 2-3 directions, which can be phrased in terms of the null complex frame 
as 

(23) 


ms e^'^'ms, 


ms e *'^m 5 , 


L 

for some function A. This will induce changes in Mjo as follows (recall that D = £ ■ d): 


Mso e*''' Msoi 

5 5 

M50 Tfgo 


(24) 

(25) 


® We will use indices i, j, ■. ■ to label both the real basis {i,j = 2,3,4) and the new basis {i,j = 4,5,5). The meaning 
should be clear from the context. 
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The last equation, in particular, implies that we can choose A to set Miq= 0 (note that the LHS 

4 

is imaginary). Moreover, eqn. NPl of Ref. |28] for ij = 45 gives M 5 o= 0 which, from the above, is 
preserved under such spins. Hence we have 

MjO= 0. (26) 

Finally, we are free to rescale t so that the geodesics with tangent (. are affinely parameterized, which 
implies Lio = 0. Note that the conditions 

Lio = 0, Ki = 0, r' = 0, Mjo= 0, (27) 

mean that we have chosen our basis to be parallelly transported along the geodesics with tangent L 


2.3 Determining a, h 


We introduce local coordinates as follows. Pick a hypersurface S transverse to I and introduce 
coordinates on S. Now assign coordinates (r, x^) to the point parameter distance r along the 
integral curve of t through the point on S with coordinates x^. Note that r is an affine parameter 
along the geodesics. We now have 

Consider eqn. NPl of Ref. [28], which in our parallelly transported basis reads 


Taking the ij = 55 component gives: 

A 

dr 

The solution is then 


D Pij — PikPkj- 


[6(1 — ia)] = — [6(1 — ia)Y 


(29) 


(30) 


6(1 — ia) = 


1 


(31) 

r -IX 

for some complex function x that does not depend on r. There is freedom in defining r in the sense 
that we can shift it by a function of the other coordinates, r ^ r + a(x^), which corresponds to 
moving the surface S used to define r. This freedom can be used to set Im(x) = 0. With y real, we 
take the real and imaginary parts of the above equation to find 

X 


a = -- 


6 = 


The r-dependence of pij is then given by 


Pij — 


0 

0 

1 


r2 _j_ ^2 


r+ix 


r-ix 

0 


In the real basis, the eigenvalues of pij are 6(1 + a^) and 6(1 ± ia). Using 
1/r and l/(r ± ix), as asserted in our Theorem. 


(32) 


(33) 


we see that these are 


2.4 The case x = 0 

If X = 0 then we have pij = (l/r)dij, so I is free of rotation and shear with non-vanishing expansion. 
This defines the Robinson-Trautman class of solutions. These solutions were studied in Ref. [20] . 
where is was proved that the only solution of this type (in 5d) is the generalized Schwarzschild 
solution (jlOp . This establishes case 3 of the Theorem. Henceforth we will assume x 7^ 0- 
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2.5 Boost weight 0 components of the Weyl tensor 

The only GHP scalars with boost weight 0 that we do not know yet are such components of the Weyl 
tensor: ^ijki, ^ij- Their r- dependence can be determined completely using the Bianchi identities, as 
we will now explain. Notice that, in five dimensions, all the information regarding the boost weight 
0 components of the Weyl tensor is encoded in for one can write 

<^ijki = -2 ^ {SikSji - 6iiSjk) (34) 


in terms of the symmetric part ‘h?- = 

The relevant Bianchi identities to determine are equations B2, B3 and B4 of Ref. [28], which 
become, in our basis, 

D^ij = - {^ik + + ^Sik) Pkj: (35) 

^^ijkl ‘^^iji^Pkl Plk) “1“ ^kiPjl ^liPjk ^kjPil ^IjPik ^ijkmPml ^ijlmPmk (36) 

and 

0 = - 2^iyPkl] + ^im[jk\Pm\l], (37) 

respectively. 

Take first the ijkl = 5455 component of (I37|l . which gives 


A ^ ix{2r + ix) ^s 

45 I o,-, ^ 45' 


r(r + 2ix) 

Taking the ijkl = 4555 component of p 6 p and substituting the expression for $45 then gives 

4 3 


(38) 


2D<^% + 1 — + 


5 

2 r 


r + ix 2r + 4ix 


^45 = 0 - 


(39) 


Comparing the previous two equations with the one obtained from the ij = 54 component of (1351) 
gives 4>45 = 0 , and hence ^>45 = <^54 = 0 . 

Now we compare the two equations that can be obtained for <^ 55 . The first comes from setting 
ij = 55 in ([35|), 

-D4>55 + 4>55 = 0. (40) 

r + ix 

The second is obtained by putting ijkl = 4545 in (l36l) . 


2T><h55 + 


5r + 3zx 

/ I ■ n ^55 - 0, 

r[r + ix) 


(41) 


One then immediately sees that consistency requires <h 55 = 0. 

Thus, the nontrivial components of can only be 4 > 44 , 4 > 55 . 
the relation 


A_ ^ 2 ixr 

^55 J.2 _j_ ^2 


$ 


44. 


Setting ijkl 


If we now substitute this into the ijkl = 4545 component of (1361) . we find 


4455 in (f37|) gives 
(42) 


Q4.2 ^2 1 

T»4>44 + , o . ^os 4»44 + -^55 = 0. 


+X^) 

On the other hand, the ij = 44 component of ()35p gives 


Zl<h44 H— (^44 T *^ 55 ) — 0- 


(43) 


(44) 
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(45) 


Comparing these we find 


and hence 


„2 _ 2 

<!>?, = " * 


55 j,2 _j_ ^2 


{r + ixf 


^55 ^2 + ^2 


<f>44. 


Substituting this back into dM]), we have 


Z)$44 = 




44 


9r 


for which the solution is simply 


<f>44 = —- 


4r 

r 2 + ;;t2 


$ 


44, 


(r2 + y2)2 ’ 

for some function [i independent of r. Therefore, the form of ^ij is 


$.. =__ 

1-] (j.2 _|_ ^2^3 


^2 _|_ ^2 
0 


0 0 
0 (r + ixY 


0 (r - ixY 


0 


(46) 


(47) 

(48) 


(49) 


One can then verify that any other component of equations (l3^ . (p6]) . ([37)1 is trivially satisfied. 

Now recall that we have chosen the normalization of ^ such that it is tangent to affinely param¬ 
eterized null geodesics. Choosing the affine parameter as one of the coordinates, as above, allows 
i to be written as = djdr. This, however, does not determine r uniquely. Previously we have 
used a shift in r by the other coordinates to make x real. Now consider the effect of a boost 1 = \l 
for some non-zero function A. The new vector field i will also be tangent to affinely parameterized 
geodesics provided that A is independent of r. Thus, we can define 


i 


A 

dr ’ 


(50) 


where f = r/A. If the analysis above were repeated using f instead of r, the optical matrix pij of £ 
would have the same form as pij if we defined x = x/^- Note that this is consistent with the fact 
that Pij transforms with boost weight -|-1, Pij = \pij- On the other hand, has boost weight 0 
and so is invariant under boosts. Therefore it would retain the same form as in eqn. (1491) by defining 
p = pjyY. But this shows that we can choose A so as to make p constant (but we can’t choose its 
sign). Dropping the hats, we can assume, without any loss of generality, that 


p = po, (51) 

where po is a constant. Note that vanishes if, and only if, pQ = 0. 


2.6 Optical matrix pC of 

Having determined the r-dependence of ^ij, the r-dependence of pC can now be completely deter¬ 
mined by using information from eqns. NP4, NP4' of Ref. |28j . In our basis, these read {A = n ■ d) 

k k ! A. 

^Pij T'llPijrT Mil PkjA Mjl Pik — PikPkj ^ij (^^) 


and ^ 

^Pij ~ ~PikPkj ~ ^ji ~ 

respectively. 


( 53 ) 
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( 54 ) 


We start by taking the ij = 44 component of (I5.sp . which gives 


A 


A'/044 = , o 2^2 “ T 

r (r2 + x2)2 4 


The solution to this equation is 


Pa4 — 


A 


Po 


Ar 


r 2 r(r^ + x^) 8 ’ 


(55) 


where A is some function independent of r. Similarly, the ij = 45 component of the same equation 
yields 


^P45 — 


to which the solution is 


P45 — 


P45 

r + ix 
B 5 


r + ix 

for some complex function that does not depend on r. Substituting this into 
one finds 

4 iB^ 

M51= — 

X 

and then the ij = 54 component of (|52ll determines 

/ _ -S 5 

Pb4 - —■ 

In addition, (I52p with ij = 55 immediately gives 

P55 = 0- 

Next, putting ij = 55 in eqn. ([53l) gives the equation 

Dn'- = ^55 , _ A 

55 r — ix (r^ + x^)^ 4 


This can be integrated to give 


^55 


Poir - *X) 1 


2(j- 2 _j_ ^2^2 r — ix 


A 5 - ^{'r + ix){r -3ix) 


(56) 

(57) 

for ij = 45, 

(58) 

(59) 

(60) 

(61) 

(62) 


where A 5 is a complex function independent of r. No further information can then be extracted 
from eqn. ([53|) . 

Consider now the ij = 44 component of (I52p . 


rf + rLii = r /?44 - 




+ 


Ar^ 


(r2+x2)2 4 ’ 


(63) 


where rf is the r-component of n“ in the coordinate basis defined by (r, x^). Similarly the ij = 55 
component gives 


rf - iAx + (r - ix)Aii = (r - ix)P55 - 2^ 2 + “ *A)^ 

T* “T X 4 

where we have used Ar = n’’. Subtracting p63p from (1641) gives the relation 

A 5 - A + iAx + ixLii + , + 5x) = 0. 

(r"' _|_ g 


(64) 


(65) 
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Taking the real part yields simply 


and hence one can write 


^55 


Re(j45) — + 

iKx IJ^o{r-ix) 1 


5Ax2 


+ 


2(^2 _j_ ^ 2^2 r — iy 


A-iF- 


Ar^ 


( 66 ) 


(67) 


where we have redefined Im(A 5 ) = —F. Thus the r-dependence of pA is now known. There is still 
some information left from eqns. (I63p and (j64|] . which will be exhausted in the next calculation. 


2.7 Determining the basis vectors 

In this section we determine the r-dependence of the basis vectors. It is sufficient to consider the 
commutators 


[£,n] = -Liii, 

[i,mi] = -Lui - pjirrij, 

together with the remaining information from eqn. (|52p and also eqn. NP3 of Ref. 
basis is (<5, = nii ■ d) 

i i 

^[j\Pi\k\ - Li\j\Pi\k\+ Mi[j\ Pi\k\+ M[kj\ Pii = 0. 

The n-components of (IHHIl reads 

Dn^ = 0 . 


Hence 


n 




where is independent of r. In order to determine n'’, we use ()63l) to write 


_ n , 

ll--y+P44-(^2+^2)2 ' 4 


Por 


Ar 
+ —. 


Substituting into (|64l) then gives immediately 


= A+ [n°(x) -T] - - 


Po 


X 8 2{r‘^ + x‘^)' 


The r-component of (1681) now gives 


_ n r _ A" - n°(x) Ar por 
11 ^ X ^ 4 (r2+x2)2' 

All the information contained in eqn. (|52p has now been used. 

Next consider the p-components of (|69p with i = 4: 

r 


The solution is simply 


k- 

m7 = 


(m^) 


0 \k 


( 68 ) 

(69) 

which in our 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 


where (m^)^ are functions independent of r. Similarly, the r-component gives the relation 

Li4 = -Dml - 


(76) 

(77) 

(78) 


II 













Substituting this into eqn. (170]) with ijk = 545 yields 

imKx) 


r 


= (r - ix) ( Dml + ^ j + y (r - zy) M 55 , 


where m^(x) ^ m^^d^x- Now setting ijk = 455 in dTOf) . gives 

4 4 

(r - ix) M 55 +{r + ix) 0 , 


(79) 


(80) 


showing that (r — ix) M 55 is purely imaginary. The last term on the RHS of eqn. (|79|) is then real. 
Taking the imaginary part of ()79p gives a simple differential equation for 


which integrates to 


d™; =-!< + 

r 


r E4 mlix) 

m^ = -^-, 

r X 


for some function 75^4 of the x^-coordinates only. Going back to (|78p . one then finds 

mlix) 


Li4 — — - 


Xr 


(81) 

(82) 

(83) 


The real part of ([7^ then determines M 55 , 


4 iE^ 

M55= 


X{r-ix)' 


(84) 


Now following the same procedure as above determines the form of m^. The //-components of (I69p 
with z = 5 imply 

Hr 


= . , 

^ r + IX 

for complex functions {m^)^ of the x^-coordinates only. Then, using the r-component to write 


(85) 


Li 5 = -Dml - 




r + IX 

and substituting into ([701) with ijk = 555 gives 

r _ 2ixml irnlix) 

"^5 ,2 + ^2 ,2 + ^2’ 

where m^x) = Integrating this equation one finds that 

_ E5{r - ix) + rm^ix) 

— -7-^—7-? 

xir + ix) 

for some complex function independent of r. Eqn. (1861) then determines L 15 : 

E5 + ml{x) 


Li5 — -- 


xir + ix) 


( 86 ) 


(87) 


( 88 ) 


(89) 


In summary, the coordinate basis components of n“, and 777,5 given by (IZ21), (ED, (ED, dSD, 
(| 8 D and ([ 88 ]) • Thus, the r-dependence of the basis vectors (and hence the metric) is fully determined. 
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i i 

2.8 Calculation of the non-GHP scalars Mji,Mjk 

The results obtained in the previous subsection regarding the r-dependence of the basis vectors 

has automatically given information about some connection components that are not GHP scalars, 

4 

namely Lu, Lu and M 55 , as given in eqns. (I7^ . (l 88 ]l . ([ 8 ^ and (l84]l . respectively. Furthermore, 

4 

our previous calculation of pT also provided M 51 , eqn. (I58|) . Here we show how to obtain the other 

i i 

components of Mji,Mjk- 

We first notice that the calculation of the previous subsection does not exhaust eqn. (l70l) . there 
are two components remaining. The first is ijk = 445, for which one finds 

M54= —, (90) 

Xr 

using previous results. The second is the ijk = 545 component, which gives 

M55= 0. (91) 

Now consider the commutator Cl of Ref. [28] applied to a GHP scalar V) of arbitrary boost 
weight b and spin weight s = 1, [t>, V). This gives two relations: a boost part (i.e. the coefficient 

of b) 

T)Lii = -$ + ^, (92) 

which is automatically satisfied, and a boost-independent part, 


D Mji= -24>A. 


(93) 


4 

Notice that this implies that M 51 is independent of r, which is consistent with the result already 

5 

known, eqn. (I58|) . On the other hand, for Mgi this equation gives 


5 


— iC 


¥oX 

(r^ + X^Y ' 


(94) 


for some unknown, real function C of the x^-coordinates only. 

One can follow the same procedure by applying the commutator C2 of [28|, [ 13 , 6 ^], to a GHP 
scalar Vi of spin weight 1. This gives the equations 


DL\i — LijPji 


and 


D Mjk= - Mji pik- 

The former is automatically satisfied, while the latter gives the additional information 

5 iD4 

MU= — 


and 


5 iD^ 

M55= -- 


r + ix 

where D 4 and ZI 5 are real and complex r-independent functions, respectively. 


(95) 

(96) 

(97) 

(98) 
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2.9 Boost weight —1 components of the Weyl tensor 

We are now in a position to take ] 3 ^- and dj-derivatives in fnll. In particnlar, we can now consider 
the boost weight 0 components of the Bianchi equation, eqns. B5, B5', B 6 and B7 of Ref. [28], and 
determine completely 

Start with eqn. B5 of [28], which reads 


- - ^j\ii) Pi\k]- 

The ijk = 545 component of this equation immediately gives that 

'^545 = 0- 

If one takes the ijk = 555 component, one finds that 

^ 2^o£^5(2r + ix)(r-zx) 4^omg(x) 

while setting ijk = 545 gives 

, ^ 2noEi{r-ix) _ 4/iom°(x) 

^545 (^2 _|_ ^2^3 ^('^2_|_^2p’ 

Substituting this result into the ijk = 455 component, one finds that 

mlix) = 0 


(99) 

( 100 ) 

( 101 ) 

( 102 ) 

(103) 


and hence 


\l/' - = — 
^545 


2fioEi{r - ix) 


(^2 ^2)3 

The only remaining independent component of (1991) is then ijk = 445, giving 


(104) 


vT;' - = _ 
^555 


2 ^ioE^r{ 2 r - ix){r - ix) 4/io(r - ix) ml{x) 


X{r‘^ + X^)"^ 


^(^2 _^^ 2)3 


Now consider eqn. B5' of 


- 26[j4>fc]i + = 2 ( ^'6uu. - T' 


ijk 


'Aj\i - Pm 


Using the previous results, the ijk = 445 component gives the condition 

"m^x) = -E5- 

Eqns. (|10ip and (|105p then imply 


4^445 — — 


2ijioE^{r - ix) 


and 


'k' - = 

^555 


{A + X^)^ 
2 i^ioE^{r - 2ix){r - ixf 


(A + x^)^ ’ 

respectively. Of course, the symmetries of determine automatically 

\I/' - = 

455 545 545 ^j.2 _j_ ^2^3 ' 


(105) 

(106) 

(107) 

(108) 

(109) 

( 110 ) 


One can then verify that all other components of (I106p are automatically satisfied, as well as eqns. B 6 
and B7 of Ref. [28] . 
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2.10 Equation NP2' 

Eqn. NP2' of Ref. [28], which reduces to 



II 

( 111 ) 

This can be integrated to determine k[. 

. The result is 


K4 

II 

+ 

to 

to 

( 112 ) 

Kg 

ipoE^ir-ixf 

(^2 + ^2)3 > 

(113) 


where G 4 (real) and G 5 (complex) are functions depending on only. 

2.11 Differential and algebraic constraints 

The r-independent integration functions appearing in the various expressions in the previous sec¬ 
tions are not completely independent. There is still information contained in the commutators 
[n,mj\ ,\mi,mj\, as well as in the commutators of GHP derivatives [^, 61 ] and applied to 

some GHP scalar Vi of spin weight 1, and this information can be used to place algebraic and 
differential constraints on the r-independent functions found above. Later on, when we introduce 
coordinates, we will be interested in the symmetries and Killing fields admitted by the solutions 
considered here. In order to study these symmetries, we will need to know the derivatives of those 
functions along all basis vectors. 

Gonsider first the commutator 


[mi^nij] = + 2p\ij^n + 2 M[ij] nik (114) 

for f = 4,j = 5. The /r-components of this equation give 

X [m 4 , mg] = -E^ml -t- {ixDi - iE^) mg. (115) 

On the other hand, the r-component can be brought to the form p{r) = 0, where p{r) is a polynomial 
in r with coefficients depending on the x^-coordinates. It is clear that the coefficient of each power 
of r must then vanish, resulting in the independent equations 

0 = x"i 4 (-E^s) - 2 x^.B5 - ixD^E^ + lE^E^, (116) 

0 = -ixml{E4,)+ iE4E5 + 3x^B5. (117) 

Since the action of m^ and mg on x is known, eqns. (jlO.SI) . (|107p . one can apply their commutator 
(|115p to X to find 

ml{E^) = iDiE^-i^^. (118) 

X 

Gomparison with (|116ll yields immediately 


R5 — 0. 


(119) 


Hence the only non-vanishing components of pL are / 044 ,/Ogg (and its complex conjugate). Using the 
above result, eqn. (|117l) reduces to 


mg {E^) 


E^E^ 

X 


( 120 ) 
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(121) 


An additional algebraic constraint can be obtained from eqn. NP3^ of [28] , 

^\j\Pi\k\ ~ ~ 2^ijk- 

By taking the ijk = 455 component and using = 0, one finds 



( 122 ) 


Now consider [ 771 , 5 , 7715 ] from eqn. (|114ll . Similarly to the calculation for [ 7774 , 7775 ], one gets the 
following independent equations: 

X [ 7775 , 7775 ] = 27x^77° — 2iE/^rrbl + ixD^m^ + ixP^sm® (123) 


and 

777 ° {E,) = ^ + 2{A + iF)x - iD,E^ - 

4 X 

If we act with [7775,717,5] on x and use (I 124 D . we find 

n°{x) = 2 P 

and hence we know all derivatives of x- 
Also, we have the commutators 

[?7,7i7j] = + Lun - ( Mil +p'ji I rrij. 


(124) 


(125) 


(126) 


From the 7 = 4 component one finds the relations 



7174 (F) = 0, 

(127) 


0 = -2FF4 + xn° (F 4 ) - X" 74 (A), 

(128) 

and 

X [n°, 777 ^] = Eml, 

(129) 

while the 7 

= 5 component gives 



0 = 777 °(F) - ^ + 77 ° (F5) - + G5X - iCE^, 

X ^ 

(130) 


0 = xml{A) 27FF5 + 7xn° (F5) + ^ iG^x^ + GE^x 

(131) 

and 

X [ 77 °, 777 ^]= tGx + E+ ^ ml- 

(132) 


Apart from (I127p . which gives directly the derivative of E along m^, the other relations determine 
only a combination of derivatives. One can, however, gain more information from other equations. 
We start by going back to eqn. NP3^, (I12ip . Taking the ijk = 445 component one finds 


07 iFEi, 

m,{A) = - - -—. 

Substituting into (| 13 ip gives 

77 ° {E5) = G,x + iCE, + ^^. 


(133) 


(134) 
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Now using this in (ll.SOp gives 


ml 


{F)={ 


( iAEk 


- 2 G 5 X- 


FE 5 

X 


Setting ijk = 545 in (|121ll and using ()127p . one finds 


ml{A) = 


FE 4 

X 


which can now be used in eqn. (jl28p to obtain 


(E^) = 


3FE4 

X 


(135) 


(136) 


(137) 


No additional information can then be obtained from eqn. NP3', (I12ip . 

Another source of information is the commutator of GHP derivatives Applying this to a 

GHP scalar V) with spin weight 1, one obtains two equations, just as in the study of Vi and 
[]3,6j] Vj carried out before. The boost part gives 


^ / A 

^lj\Fi\i] = -FiiP[ij] - Lik M[ij] +Pk[i\Pk\j] + 
which is automatically satisfied using previous results. The boost-independent part is 


(138) 


jl JXIjk 


— PikP'jl + PilP'jk + PjkPil — PjlPik + 2 Mjl P[kl] 
i P i P ip A 

T MpkMjl MplMjk T2 jpM^kl] ^ijkl {^ik^jl ^il^jk) • (139) 


Setting ijkl = 4545 gives 
while the ijkl = 4545 gives 


"I5 (^5) = -iD^E^, 

ml {Ei) = -3Fx- 


(140) 

(141) 


There is now enough information to apply the commutators [n^, 

considered here. Before we used [m^jTTig] (x) to find eqn. (I125p . 
relation can then be obtained by considering (x), namely 


m. 


01 


l J 5 

Anot’ 


ml, m^ 


to the functions 
rer non-trivial, algebraic 


SiAE^x 

4 


^G'sX — iCE^ — 0. 


(142) 


Applying the same ideas to A, one finds that (A) gives 


E 4 




(143) 


It turns out that the term in brackets always vanishes, for suppose that E 4 = 0. Then (I14ip implies 
that F = 0, which in turn, from (|135p . implies 


G5 = 


iAE^ 

4 


(144) 


Thus, the latter is true irrespective of whether E 4 = D ox E 4 ^ 0. Using this in (I142h . one finds 
another important algebraic constraint. 





= 0 . 


(145) 
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In turn, (^4) gives 


3AFx 2AF 


n^(A) = + 


X 


Finally, (A) implies either E 4 = 0 ov 


n°(F) = 


2 F^ 

X 


(146) 


(147) 


However, as just discussed above, £^4 = 0 implies £ = 0, which is consistent with (| 14711 . Therefore, 
one can safely take (I147p to be always true. 


2.12 Equation NPl' 

Now that we have the information about all derivatives of the functions involved, we can finally 
consider eqn. NPl' of Ref. [28j to calculate Ifb, the Weyl tensor components with boost weight —2. 
In our basis, eqn. NPl' of |28j reads 


k k / k 

^Pij ~ + 2K[Lij + k'^ Mij —p'kj Mil —p'ik \ p'kj+ Mji 




(148) 


Taking the trace of this equation and recalling that Pb = Q, one finds 

£ = 0 . 


(149) 


Then, the independent components of Pb can be obtained by putting ij = 44,45,55 in eqn. (|148ll . 
giving 




44 




45 


O', 


55 


IfiO (-E| — E5E5) 
(r^ + 

GippE^E^jr - ixf 

(^2 _j_ ^2j4 

6/io£g(r-ix)^ 

(’^2 + y 2)5 


(150) 

(151) 

(152) 


with Hgg = —H 44/2 following from the traceless condition. No further information can be obtained 
from eqn. (jl48p . 

Note that all Weyl curvature components (i.e. ^ij, Tb^ and Pb) are proportional to pp. Fur¬ 
thermore, 7 ^ 0 if /iQ 7 ^ 0. Hence the solution is conformally flat if, and only if, ^0 = 0- 


2.13 Further algebraic constraints 

Notice that, with £ = 0, we have 


n°(x) = 0 , m^(x) = 0 , mg(x) = -£ 5 , 


n^(A) = 0, m 4 (A) = 0, “ 


n°(£4) = 0, rnfl{E4) = 0, mg(£4) = 


A£5X 

4 ’ 

X 


With this information, we find 


£5 


£5 


VaX =-— ("l5)a - , • 

r — IX r + IX 




(153) 


(154) 
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When computing one then hnds 


V 


a 



= 0 , 


showing that the term in brackets must be constant, 


A = Ao + 



for some constant Aq. Similarly, the equation for Va £'4 becomes 


(155) 


(156) 


VaixEi) = 0 (157) 

and hence 

Ei = — (158) 

X 

for some constant Eq. 

We have now finished calculating all the components of the connection, including non-GHP 
scalars, and the components of the Weyl tensor. Furthermore, we have found the r-dependence 
of the basis vectors and several algebraic and differential constraints involving the r-independent 
functions. 


2.14 Null rotations and type D 

Consider a null rotation (fTTll about £. Choosing 

Eo iE^ 

ZA = -, ^5 =-, 

Xr r + ix 

one finds that the following CHP scalars transform according to 

k' !-)• 0, !-)• 0, nP !-)• 0 

and 

/ 1 + {a'f 0 0 \ 

I 0 0 1 - ia' , 

\ 0 1 + ia' 0 / 

where 

/ X 
a = -, 

r 

U = _ Ar(r^ - y^) E^ 

2(^2 _|_ ^ 2)2 ^2 _|_ ^2 8 ( j ..2 _j_ ^ 2 ^ 2 rx ^{ r ‘^+ 

We have thus found a basis in which only the boost weight 0 components of the Weyl tensor are 
non-zero, and where both £ and n are geodesic with corresponding optical matrices in their canonical 
form. This shows that all solutions considered here are of type D. 

If 6' 7 ^ 0 then the optical matrix of the second multiple WAND, has rank 3 (in agreement with 
an argument of Ref. mi)- If 6' = 0 then the optical matrix of the second multiple WAND vanishes 
identically in which case the solution belongs to the Kundt family. The condition for this is 

b' = 0 ^ e^E^ = ^, ylo = F;o = A = 0. (164) 

We will study this special case in more detail in section |H 


(159) 


(160) 


(161) 


+ 


- X^) + X^) 


E^E^r 

2 I .,,2'i2 • 


(162) 

(163) 
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3 Coordinate basis calculations 


3.1 Integrable submanifolds 

So far we have kept the coordinates arbitrary. We will now show how the results derived above 
lead to a canonical way of choosing these coordinates. From 


L L L 

Mab= -KiLh - PijL{mj)b - Pijna{mj)b+ Mji {mj)Jb+ Mjk {mj)a{mk)b, (165) 
one finds, using i ?5 = 0 , that 


dms = — 


5 4 5 5 

(/’ 55 + + (ikfss + M54)m4+ M55 


A ms. 


(166) 


Thus, ms A dms = 0. This is the integrability condition for the existence of a complex function z 
such that ms = m® oc dz. Thus, we write 


ms = A4dz 


(167) 


for some function A4. 

The distribution {£,n,m 4 } is tangent to surfaces of constant z and hence integrable. This can 
also be seen from the commutators (| 68 |) . (f69]l . (11261) or, more explicitly, 

(168) 

(169) 

(170) 


[£,n] = -Lni, 

[£, m4] = —Li4i — P44m4, 

[n, m4] = —K 47 + Li4n — p'44m4. 

We also have the commutator (jl29p (using F = 0) 

[n°,m4] = 0. 

Thus { 11 ^, 171 ^} is integrable and we can choose coordinates {r,u,x,z,z) so that 


n° = - 


du 


d 


d 


mA = N I ^ - 


dx 


' du 


(171) 


(172) 


for real functions N 7 ^ 0, Lx independent of r. From (ms)^ = 0 we have m^ = 0 and hence 


ms = M 


A 

dz 


-Y 




for complex functions M 7 ^ 0, Y, independent of r. We now have 

n_ d _ rjd _ d 

^ dr' ^ ^ dr du 


where 


H = rF = Aq- 4(r-^ - X^) - 


Po 


(173) 


(174) 


(175) 


2(r2 + x^) 

depends on r. Now define I 4 (real) and Vz (complex), both independent of r,hy E 4 = Eq/x = —NVx 
and E^ = —iM{Vz — YVx). We then have 


NVx d N ( d d 

m4 = -H-^-Ta,— 

r or r \ ox du 


(176) 
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(177) 


M(14 - YV^) d M 

"^5 =- 

r + ^x or r + IX 
From the inner products we then find 


A 

dz 


-Y 


A 

dx 


-L. 


A 

' du 


- L 


A' 

' du 


i 

n 

7714 


ms 


—du — Lxdx — Lzdz — L^dz, 
dr — Hi + 14dx + Vzdz + Vzdz, 
^ (dx + Ydz + Y dz) , 


from which it is easy to write down the metric. 

Recall that n^{x) = nr^ix) = 0. These imply that 


duX = dxX = 0 , 


hence y = xiz,z). Now, im}l{x) 


—E^ reduces to 


dzX = 


M’ 


while the commutator = 0 gives 


duN = duLx = 0 . 


(178) 

(179) 

(180) 

(181) 


(182) 


(183) 


(184) 


From (11581) we then have duVx = 0. The result above implies that Lx = Lx{x, z, z). Note that we 
have the residual coordinate freedom u ^ u' = u + h{x, z, z), which has the effect 

Lx L'^ = Lx — dxh, Lz —>■ = Lz — dzh. (135) 

We can therefore choose h appropriately to set L'^ = 0, and then drop the primes. Henceforth, 



II 

o 


(186) 

The expression 

[n\ml]=i{c + ^'^ 

mg 

(187) 

reduces to 


C + M. 


duLz 

= duY = 0, duM = -i ( 

(188) 

The latter equation implies 

= 0, 


(189) 

and 

^ Ex ^ i \'nP{M) n 

^“4+2 M 

0(M)' 

M 

(190) 


where we used the fact that duM/M is purely imaginary to write (7 in a form that is manifestly 
real. If we now recall that 'nP{E^) = 0, we find 


0 = duE^ = -iMduVz, (191) 

where we make use of (jl45p . Hence, duVz = 0. Thus, we conclude that all functions appearing in 
the metric 

9ab — ^art-b T T (rUi)a{rUi)b (192) 
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are independent of u, and hence d/du is a Killing vector field. Notice, in particular, that the metric 
depends on M only through |Mp. 

The commutator 


[ml,m^] = 


E 5 


- niA + i , „ 

X V 


£>4- 


Er 


ml 


(193) 


gives 


dxEz — 
dzixN) = 

ml{M) = 

The latter equation can be used to show that 


dzLx = 0 , 

(194) 

xYdxN - xNdxY, 

(195) 

iM (^Da - . 

(196) 


(197) 


5.|M|2 = 0, 


as well as to write ZI 4 as 


ZI 4 


NVa, , i \ml{M) ml{M) 
^ 2 |_ M M 


(198) 


Moreover, there is still freedom in redefining x by transforming x —>■ x'{x,z,z). This changes N to 
N' = Ndx'/dx. Since N = N{x, z, z), we can use this to impose the condition N' = 1/x- Dropping 
the primes, N = 1/x, which implies that the LHS of (11951) vanishes, as well as dxN. Thus (11951) 
simply reduces to 

dxY = 0. (199) 


Notice also that K 4 = Eq/x = —NVx implies that I 4 = —Eq = constant. The equation for 


ml{E^) = W^E^ -( 200 ) 

then reduces to 

dxVz = 0 . ( 201 ) 

Therefore, similarly to the argument above for d/du, we have shown that d/dx is also a Killing 
vector held for these solutions. The function M depends on u and x only through a phase, which 
may be eliminated with an r-independent spin transformation of the form (|23p . 

Consider now the expression for [ml, , 



[ml, mg] - 2 zxn° "^4 + ^D^^l + iD^ml, 

( 202 ) 

which reduces to 

2 ix = \M\^dzLz-dzL,), 

(203) 


- au’) . 

(204) 


n iml{M) 

= M ■ 

(205) 

With the above equation for D 5 , the equation 



777-5 (Ks) = -iD^E^ 

(206) 

can be written as 

ml{E^M) = 0 . 

(207) 
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But if we recall that = —MdzX^ see (|183j) . the quantity in brackets is —\M\'^dzX, which is a 
function of z, z only. The equation above then shows that it must actually be a function of z only: 

E 5 M = -\M\^dzX = k, (208) 

for some analytic function k = k{z). If we now consider a holomorphic transformation 2 ; ^ 
the effect is to change M —)■ M' = Mdz'/dz and hence k ^ k' = kdz'/dz. We can therefore use this 
transformation to set k' to —1 or 0. Dropping the primes, we have shown that we can write 

E^M = —\M\'^dzX = ko = constant G {—1,0}. (209) 

We now see that dzX = dzX- If we now write z = z\ + iz 2 -, we see that 

d2X^^ = 0 , ( 210 ) 

hence y = xi^i)- We can now distinguish between two different cases. 


3.2 Case 1: dx ^ 0 

We will now prove case 1 of our Theorem, dx ^ 0 implies dzX 7 ^ 0 so we must have feo = —1. The 
expression 

Ay^ - E"^ 

^5(-®5) = 2 A 0 X H— - iD^E^ -1, (211) 

which is derived from eqn. (I124p and previous results, reduces to 

Using dzX = 5ix/2 = l/|Afp, we can rewrite this as 


Since the quantity in brackets is a function of zi only, it follows that it must be a constant Cq. Thus, 


dx 

dzi 


Co - ^ - 2Aox' - ^ 
X^ 4 


= ^’(x). 


(214) 


We can therefore use x ^.s a coordinate rather than zi, with the transformation rule given by (I214p . 
Note that dix = 2/|Mp implies 



(215) 


from which it follows that x must lie in a range for which P{x) > 0 . 

The only quantities that we do not yet know are the functions Lz and Y, which are determined 
by eqns. (I203p and ([204]), respectively. Both have the same structure: 


\MmdzE-dzE)=gix), 


(216) 


where E denotes Lz or Y and the RHS is a known function. Notice that 

d{Fdz + Fdz) = {dzE - dzE) dz A dz 


(217) 


and the quantity in brackets on the LHS is precisely the combination in which both Lz and Y appear 
in the metric. If a particular solution Edz + Edz is found, any other solution will differ from this 
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by a gradient da, where a = a(z,z). In the case of L^, this can be absorbed by defining a new 
coordinate u' = u — a{z,z), which does not change any other quantity and, in particular, = 0 is 
maintained. Similarly, defining x' = x — a(z, z) eliminates the gradient da from the expression for 
Y (notice that dx'jdx = 1, and hence every other quantity is unchanged). Hence all we need to do 
is to find particular solutions for and Y. 

Consider first the equation for L^. If we search for a particular solution in which depends 
only on x, (I 2 n.‘lh becomes 

dhz dLz 


This can be solved, in particular, for 


dx dx 


Lz = - 


= 2ix- 




Similarly, if we look for a solution for Y such that T = Y{x), (I204p reduces to 

‘HEq 


A solution is 


dY dY 

dx dx 

Y = 


X^ 


iEp 

2X^' 


Then, from the definition E^ = —iM{Vz — YVx) and from (|209ll . (|2I4I) . one finds 


—K-f 




(218) 

(219) 

( 220 ) 
( 221 ) 

( 222 ) 


Thus, all functions appearing in the metric ()192l) depend only on x (equivalently z\) and r. In 
addition to d/du and d/dx, we now find that djdz^ is also a Killing vector field for this metric. 
Using our definitions and results above, the metric becomes 


d^ = —2(dri + x^d.Z2) 


dr + ff(du + x^dz 2 ) - Ep I dx - ^dz 2 ) + Edz 2 


X" 


+r^X^ - ^dz2'j + (r^ + X^) 

Making the definition y = Z 2 this gives the metric ([H in the statement of our Theorem. 


(223) 


3.3 Case 2: dx = 0 

In this case, x is constant. Eqn. (12091) gives E^ = kp = 0. The equation for m^{E^) is then trivial, 
but the one for 7 fi^{E^), eqn. (j 211 h . becomes 


A — -^0 
Ap — 


El A 


2x4 4^ ’ 


(224) 


which gives an algebraic constraint involving these constants. There is a non-trivial component of 
(11391) that has not yet been considered. The ijkl = 5555 component gives 

- - 
0 = i~x^ [m^ ( 1 ) 5 ) - mg (Dg)] + 2AoX^ + ^X^ - 2C'x^ + ^ + 2EoT>4 + 2D^D^x^. (225) 

X 

In general, using our expressions for C, 1 ) 4 , this gives a second-order differential equation involv¬ 
ing M, M which can be put in the form (using (j224D ) 


= ^ + 2Ax^ 


(226) 
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where is the Ricci scalar of the two-dimensional metric 

= pIp''"'’"- 

Hence this two-dimensional metric has constant curvature. A volume form for this two-dimensional 
metric is 

^ (228) 

Eqns. (j203ll and (j204[) then become 

(i{Lzdz + L^dz) = 2xe^‘^\ (229) 

- 2 

d{Ydz + Ydz) = (230) 

X 

respectively. If we define a one-form A = Az{z, z)dz + Az{z, z)dz by 

e(2) = d.A, (231) 


then particular solutions to the equations above are 


Lzdz + Lzdz = 2xA, 


Ydz + Ydz 



(232) 


As in case 1, any other solution would differ from these by some gradients da{z,z) and d/3(z,z), 
respectively. These can be absorbed into u and x using the residual coordinate freedom u' = 
u + a(z, z), x' = X + /3(z, z), which preserve all quantities fixed above. Then, from £'5 = —iM{Vz — 
YVx) = 0 one finds I 4 = —EqY, and hence 


The metric is then 
ds^ 


Vzdz + Vzdz 


2_^ 

■v3 


A. 


—2{du + 2xA) 


dr + H {du -|- 2xA) — £0 ( dx -|- 


+r‘^X^ ^dx • 


2Eq 


A 


(233) 


(234) 


We now transform to arbitrary real coordinates y°‘{z,z) on the 2d part of the metric, so that 
g(2) = ha/ 3 {y)dy°‘dy^, and this gives the metric ([7|) in the statement of our Theorem. 

The metric (j234l) has a scaling symmetry analogous to (I12p . For A 7 ^ 0 we can perform the 
coordinate transformation 


^ \ / 

u = — r = Xr 

A 


X = 


A2 


(235) 


and the metric takes the same form as before but now with the constants rescaled as 


/ X 
^ =A 


/ _ 


p/ _ Eq 


(236) 


and ( 7 ( 2 ) replaced by 


A 2 g(( 2 ) and A replaced by A! = A2aI. 


4 Relation to Kerr-de Sitter 

In this section we will perform coordinate transformations to demonstrate how the metrics (jH) and 
(0 are related to the 5d Kerr-de Sitter solution. 
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4.1 Case 1: dx ^ 0 

If we define the 1-forms 


= du + x^dy, cr^ = dx-rdy, = dy 




then the metric (jl|) can be written 


2 2 

ds^ = —2a^dr -|- dx^ + hijcr^a^, 


where 


P 


-2H Eo -P 

„2,,2 


hij = Eo Px 0 


(237) 


(238) 


(239) 


-P 0 {r^ + x^)P 

If we let = {u, X, y} then the metric can be written 


ds^ = —2aidx^dr + gjjdx^dx'^ + 




P 


-dx 


where 

0,1 — fjj, gjj — hijOjCTj. 
Now consider a change of coordinates 

dx^ = dy^ -|- ^^(r)dr 


(240) 

(241) 

(242) 


for some functions A^{r). We want to choose the functions to eliminate dy^dr terms from the 
resulting metric. This requires 

A^ = g^^aj. (243) 

We can only do this if the RHS above is independent of y. We find thall§ 

,r .2 


A^ = - 


F{ry 


= 


Eo 


r'^E{r) ’ 


713 = - 


1 


E{r) 


(244) 


where 


E{r) = (r^ + x^) ( 

A 




p2 

2H+ ^0 


+ 


P 


^ 2^2 ^2 ^2 

= —y-r‘^ + ‘^Aor‘^ — iJ,o + Co + ^^ 


r2 ■ 


(245) 


So the RHS of (I243jl is indeed independent of x and the coordinate transformation is permissible 
provided we work in a region where F{r) y 0. Note that 


F{r) = P{ir) - /xq- 


(246) 


The metric in the new coordinates is 


ds^ 


{r^ + X^) 


dj .2 ^ ^^2 

Fyy^Pix)_ 


+ hijPuy 


(247) 


®To do this computation it is convenient to define yi = d/du, y2 = d/dx and 773 = 9/922 — x^dfdu -|- (Eq j'>d)djdx 
so that cr*(yj) = We then have where h*-’ is the inverse of hij. 
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where is defined by replacing with 


in (T*. Now using (j245p to eliminate H we can write 


+ -P^ [„1 - (r= + + -P [Eo-'' + rW]\ 


" ^2 _j_ ^2 

For Eq ^ 0 we define 

and the metric is 

ds^ = {r^+x^) 
+ 


P + 




y^ = V’o, = Eo'ip2, / = ^i, 

F{r 


d ^2 ^ ^^2 

[FiP^Pi^ 


(dz/'o + X^dV^i)^ 


r 2 + y 2 


[^^0 “ r^dz/^i]^ + [dV'o + (X^ - r^)dV’i + r^x^dV’ 2 ]' 
^ _j_ X ^ X 


(248) 

(249) 


(250) 


If we now define xi = x and X 2 = ir then this is the Kerr-de Sitter solution [9] with two non-zero, 
unequal, spin parameters, as written in eqn. (22) of Ref. |12jFI 

The case Sq = 0 corresponds to the Kerr-de Sitter metric with a single non-vanishing spin 
parameter. It can be obtained from the solution as written in Ref. m by defining 'tp 2 = 1 P 2 /Eq and 
taking the limit Kq —>■ 0 . 

Now we return to the special case F{r) = 0, for which the above coordinate transformation no 
longer works. This condition can be understood geometrically as follows. The metric admits a 3d 
isometry group (associated to the Killing fields d/du, djdx, djdy) with 3d orbits. The condition 
F{r) = 0 is the condition for these orbits to be null everywhereo Note that 


F{r) = 0 A = ^0 = Ko = 0 Cq = Ao- 


(251) 


This is precisely the condition (I164p for the second multiple WAND to have vanishing expansion, 
rotation and shear, i.e., for the spacetime to be KundtH 
In this special case, the metric simplifies to 

2 2 

ds^ = —2dr idu -|- y^dy) H— „ „ idu — r^dy)^ -|- r‘^x^dx‘^ + ^ ^ dx^, (252) 

^ r2 -H X yo 

where //q > 0 follows from T’(x) > 0- (The scaling symmetry (fT^ could be used to set yo = 1.) For 
this metric, the second multiple WAND is 


k = 


1 

P + 


■9 91 

P - 1 - 

du dy\ ’ 


(253) 


which obeys ka = —{dr)a so surfaces of constant r are all null. A surface r = ro is a Killing horizon 
of the vector field r^d/du + djdy. The surface gravity vanishes, so this spacetime is foliated by 
degenerate Killing horizons. The metric is smooth at x = 0 (for r 7 ^ 0) provided that x is identified 
with period 2Trjx/JIo, which implies that cross-sections of the Killing horizons have topology 
(assuming the coordinates u,y are non-compact). The curvature diverges at r = x = 0 (see (H9]) l. 
The curvature vanishes as r — 00 so the spacetime is asymptotically locally flat. 

Obviously this special case is a limit of the generic case with F{r) 7 ^ 0. One can obtain the 
solution as a limit of a regular Myers-Perry black hole solution. Start from a single spinning Myers- 
Perry black hole, which has A = Kq = 0 and Aq > 0 and Cq > 0. Using the scaling freedom (fT^ 

^ The parameters of Ref. |12| are given in terms of our parameters by ci = 2To, C 2 = A/4, c = Eq, bi = ( 70 / 2 , 
62 = Co/2 — /io/2. 

®This is analogous to the 4d solution given in eqn. (24.21) of [I], which can be obtained as a limit of the Kerr-NUT 
solution. 

®To see this, note that (I2U9II implies \Ep = 1/|M|^ and we have l/|Mp = P(x)/2 = Co/2 using Ao = Eo — 0. 
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we set Aq = 1/2. Then the MP spin parameter is a = '/Co and the mass parameter is /io- It 
has a regular horizon provided no > c/- The extremal solution with ^o = o? is nakedly singular. 
To take the limit, perform the rescaling (I12|) and take the limit A —?• oo with /tq/A^ fixed and 
(^0 ~ ^ 0- This corresponds to scaling the Myers-Perry black hole towards the extremal 

solution whilst simultaneously taking its mass to infinity!^ 

4.2 Case 2: dx = 0 

In the metric ([7|) written as in (I234p . set 


and let 

The metric becomes 

= —2cjdr — G{r)a‘^ + / 


x' Eo 

X = - 1 - -rU 

X x^ 


a = du + 2xA. 


dx' + ^ (l + ^ \ a 


+ + X^)9^‘^\ 


(254) 

(255) 

(256) 


where 


G{r) = 2H{r) + 


El 


El , X^\ 2 , 2 ^ do 

~ “(^“1 2 ) ~ Tv +x) 2T—2' 

V rV 4 A + x^ 

Now let 

du = dp + A{r)dr, dx^ = d(? + B{r)dr 

and choose A{r), B{r) to eliminate dpdr and dgdr terms from the metric. This giveJ^ 


A{r) = — 


G/y 


B{r) = 1 + 


A 


En 


J x^G{r) 


(257) 

(258) 

(259) 


and the metric becomes 


ds^ = 


d ^2 

W) 


- G/y + / 


dq + 


Eq 

X^ 


1 + 


X 


+ y+xyg‘'^\ 


(260) 


where 

u = dp + 2xA. (261) 

The local symmetries of , which has constant curvature, extend to symmetries of the full metric 
(j26Up . hence this metric is cohomogeneity-1, where the surfaces of (local) homogeneity are surfaces 
of constant r. 

If Pio = 0 then the Killing vector field d/dq is hypersurface orthogonal. The scaling symmetry 
(j236p can be used to eliminate one parameter, so these form a 1-parameter family. These solutions 
were discussed in Ref. [29]. 

^°The limit involves “zooming in” on the equatorial plane of the black hole, inside the ergoregion, which explains 
why the limiting metric is non-stationary. 

Here we assume that we are working in a region with G{r) yf 0. If G(r) = 0 then po = 0 in which case the Weyl 
tensor vanishes so the solution is conformally fiat. 


28 

















To analyse the metric (j260p . consider the case for which > 0. Perform the coordinate 
transformation 

P = X^-^q, (262) 

where the coefficient of the second term is chosen to eliminate the O(r^) terms in the component 
of the resulting metric. After defining new coordinates p = + and t = XQ, the metric is 


^ ^ (dV' + 2A- n{p)dtf + p‘^g^‘^\ 


(263) 


where 


and 


Pip) 



" 


_ M 

p 2 


+ 


POX 



nip) 


^PqEq 

p^7(p)i?(2) 


(264) 

(265) 

(266) 


We can use the scaling freedom (|236l) to set = g go Q implies Ax2/4 < 1. If this inequality is 
strict (which is always the case for A < 0) then the above metric is the Kerr-de Sitter solution with 
equal rotation parameters as written in m, where the 2 parameters (spin and mass) are 


a = 


X 




(267) 


If the inequality is saturated, i.e., Ax^/4 = 1 (only possible for A > 0) then we have Eq = 0. We see 
that this solution can be obtained by taking a limit of the Kerr-de Sitter solution in which a ^ oo, 
M —)• 0 with Ma^ approaching a finite limit. 

If A > 0 then > Q is the only possibility. If A = 0 then it is also possible to have 7?^^^ = 0, 
which requires Eq = 0. If A < 0 then we can have = q (with Eq / 0) or < 0. Note that the 
coordinate transformation (I262p that brings the metric to the form (I263p is valid also if 7?^^^ < 0. 
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